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Abstract. The high-energy Universe is potentially a great laboratory for searching new light 
bosons such as axion-like particles (ALPs). Cosmic sources are indeed the scene of violent 
phenomena that involve strong magnetic field and/or very long baselines, where the effects of 
i i the mixing of photons with ALPs could lead to observable effects. Two examples are archetypal 

Py of this fact, that are the Universe opacity to gamma-rays and the imprints of astrophysical 

magnetic turbulence in the energy spectra of high-energy sources. In the first case, hints for the 
existence of ALPs can be proposed whereas the second one is used to put constraints on the 
fQ ALP mass and coupling to photons. 
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, ^, 1. Motivations for axion-like particle searches 



The Standard Model of particle physics reproduces data incredibly well pQ. Some of its 



foundations are however not completely understood, like for example the absence of CP violation 
in quantum chromodynamics (QCD). The most general QCD lagrangian includes a complex 
phase term which -if not exactly zero- induces CP violation. The non-observation of even a 
very small electric dipole moment for the neutron [2j implies that this phase is smaller than 
10 -11 . This fact looks unnatural an calls for a explanation. A possible one is given by making 
this phase a dynamical field which value is driven to zero by the action of its classical potential. 
This is made possible by the introduction of a new U(l) global symmetry which is spontaneously 
broken at some scale / (this is the so-called Peccei-Quinn symmetry [3J). A new particle that 
is called the axion is then predicted as an associated pseudo Nambu-Goldstone boson [4, 5j. 
In the original idea of Peccei and Quinn, / was of the order of the electroweak scale (EW), 
inducing a mass of ^100 keV for the axion, which was quickly ruled out (see [6J [7] for details). 
Then it was proposed that / was much greater than the EW scale, leading to a very light 
and weakly interacting axion (dubbed "invisible axion"). Axions are experimentally searched 
through their coupling to photons, from the Sun [HI |9], assuming they make up the galactic 
dark matter [10J or with LASERs [TT] . Beyond the case of the strong CP problem and axions, 
axion-like particles (ALPs) appear in many models of physics beyond the Standard Model such 
as string theory |T2l El E] as pseudo Nambu-Goldstone bosons associated to the breaking of 
U(l) symmetries. The properties of the associated particles are similar to that of axions, but in 
general their mass and coupling to photons are not related, making the corresponding parameter 



space larger. For the general ALP case, the interaction term with photons is 



C = --g ia F^F^a = g ia E ■ B a 



(1) 



where g ia is the dimensionful coupling between photons and ALPs, F is the electromagnetic 
tensor and a is the ALP field. On the right hand-side of Eq. [l| the FF term is expressed as a 
scalar product of the photon electric field and the magnetic field, revealing the fact that ALPs 
can couple to photons in the presence of an external magnetic field. 

In the present article, the use of natural (astrophysical) environments to search for ALPs 
is emphasized, like during the propagation of very high energy gamma-rays over cosmological 
distances, and the effect of astrophysical magnetic turbulence on high-energy photon source 
spectra. First the conventional view of the problem of the opacity of the Universe to gamma- 
rays is presented, with the discussion of a possible indication for an anomalously transparent 
Universe. Although the possible tensions can be solved in a conventional way, they can also be 
released by invoking ALPs mixing with photons. Then it is shown that this observable could 
be used as a signature when one tries to make a discovery, but that some uncertainties prevent 
from using it to derive robust constraints. It is then shown that constraints can be obtained by 
considering the effect of magnetic turbulence around the sources and finally some examples of 
constraints are given. 

2. The transparency of the Universe and ALPs 

2.1. The conventional view of the Universe opacity to gamma-rays 

Very high energy photons (with ^TeV energies) traveling through the intergalactic medium 
encounter different populations of background radiations. The most numerous type of 
background photons belong to the Cosmic Microwave Background (CMB) and a second 
population is the extragalactic background light (EBL). The latter has a double bump structure, 
that comes from direct starlight and emission re-processed by interstellar dust in the infrared 
band as sketched on Fig. [Tj Direct measurement of the EBL is very difficult because of 
foregrounds and infrared radiation by the instruments. 
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Figure 1. Left: Spectral energy density of the cosmic background photons including the CMB 
and the EBL (inspired from [15J, with permission of the authors). Right: Illustration of the pair 
production process. 



TeV gamma-ray astronomy is sensitive to the EBL density and spectrum as it is responsible 
for the attenuation of extragalactic source fluxes at high energy. The reason for that is the 
pair production process 7TeV7EBL —> e + e~, for which the threshold lies at TeV energies in 
the terrestrial frame. For instance considering £"ebl ^0.1 eV, the threshold energy satisfying 
^th^EBL > Trie (where m e is the mass of the electron) yields E^ ~ 2.6 TeV. In Fig. [I] the 
typical range of the TeV absorption range is indicated by the horizontal arrow. Because of the 
pair production process, the highest energy photons have a larger optical depth. Before 2006 
it was commonly admitted that is was very unlikely to detect TeV photons from sources above 
z ~ 0.2. The situation changed after HESS observations of two active galactic nuclei (AGNs) at 
z = 0.186 and z = 0.165. As reported in [16], when unfolded from the EBL effect, the intrinsic 
spectra of these sources were found too hard and in tension with the source models. This 
was the first indication for a Universe slightly more transparent than expected at high energy. 
Later, AGNs were observed at redshifts as high as 0.536 by MAGIC [17] and possibly 0.61 by 
HESS [18J. The Universe is indeed more transparent to gamma-rays that expected. That puzzle 
has conventional solutions, it could be for instance that spectra are actually harder, this can 
be realized for instance including hadronic components in the AGN jets or in relativistic shock 
acceleration models. The tension can be removed as well with a revision of the EBL models, 
mainly with a lower density as for instance in the model of |19| which is compatible with all 
observations. TeV observations are now even used to provide not only upper limits on the EBL 
density but actual measurements [20]. The current situation is illustrated in Fig. [2j extracted 
from [21, 22J. It represents the energy above which the absorption becomes significant (defined 
by a optical depth r(E) = 1) as a function of the redshift of the source. The lines correspond to 
models or lower limits for the EBL density and redshift evolution from JT9J [23l [23] • Constraints 
from the spectral indices of different sources are shown as arrows, and the HESS measurement 
corresponds to the blue band. One source seems to be in tension with the measurement. However 
the methods that lead to the constraint and the measurements are different as the constraints 
rely on spectral slope measurements and the measurement comes from the observation of features 
in the spectra that can be related to the EBL spectral density. A more unified approach might 
be necessary to get a definite answer on how strong the tension is. Note also that the Fermi 
collaboration did the same measurement at higher redshifts and found a good compatibility with 
EBL models [25]. 
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Figure 2. Energies corresponding to r = 
1 for different EBL models, constraints from 
very high energy gamma-ray astronomy and 1-a 
measurements from HESS (figure from [21, 22J). 
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Some studies however still claim for an anomaly, even with the lower EBL limits from [24j. 



It is the case in [26J , where the authors claim for an anomaly, with the caveat that their claim 
requires to leave out some error bars. 

2.2. How ALPs enter the game 

The lack of opacity of the Universe to gamma-rays gave rise to the idea that ALPs could be 
responsible for this effect. The basic idea is that if mixing between ALPs and photons occur, 
the beam could travel in the form of ALPs on a significant fraction of way, not producing pairs, 
as sketched on Fig. [3| If ALPs are converted back to photons before observations, this could 
lead to a more transparent Universe. 
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Figure 3. Illustration of the photon/ ALP 
oscillations in a magnetic field 



Figure 4. Illustration of the modeling of the 
extragalactic magnetic field 



To get an idea of the relevant masses and couplings for the ALPs that are invoked here, 
let us consider the most simple formalism for describing the photon/ ALP mixing. The system 
propagation is described by a Schrodinger-like equation: 



(E - id z - M) 
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where A# = g ja B t /2 describes the photon/ ALP coupling (B t is the transverse projection of the 
magnetic field), r is the optical depth related to EBL absorption and A a = —m 2 a /2E accounts 
for the ALP mass. The fact that the imaginary coefficient only applies to the photon part of 
the wavefunction leads to the change in the overall transparency. In the case of no absorption, 
the mixing matrix is diagonalized with a rotation angle 9 such that tan 29 — — 2A j g/A a . The 
resolution of the propagation equation in the propagation state basis leads to the probability of 
transition 
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The overall 1/2 coefficient in P 7 ^ a accounts for the two polarizations of the photon. A critical 
energy E c appears, that defines the energy scale at which strong mixing occurs. From the 
expression of E C1 with cosmological magnetic fields B ~ 1 nG, an ALP mass m a ~ neV and a 
coupling g ia ~ 10 -11 GeV -1 , the critical energy lies at the TeV scale. It follows that the type 



of ALPs that are concerned by the so-called transparency hint will fall in a region of low masses 
and with couplings larger than those of the corresponding axions. 

The full treatment of the transparency problem in the presence of ALPs requires a 3 x 3 
mixing matrix to account for the two polarization states for the photon, and a description of 
the magnetic field on the path from the source to the observatory. The extragalactic magnetic 
field is usually described as a patches of coherent domains of 1 Mpc size. The magnetic field 
strength is the same in all domains but from one domain to the next its orientation changes in 
a random way (see the sketch of Fig. [4]). It can be shown (see [27]) that for random orientations 
and a large number N of domains, the transition probability is reduced to 



P^a = -(l-exp(-37VPo)) 



(4) 



where Pq is the transition probability in one domain. From this expression one would expect 
to have a 1/3 drop in the energy spectrum above E c in the limit NPq ^> 1, and a flux that is 
boosted at high energy (typically above the pair-production related cutoff) as described in [28] . 
At least two facts lead to revise the above statements. First, in practice the limit NPo » 1 
is hardly realized. Second, due to the unknown nature of the magnetic field configuration, the 
prediction on the transmission has an intrinsic variance. Indeed it can happen that the ALPs 
do not convert back into photons before reaching the Earth, leading in that case to an even 
more opaque Universe. This is nicely illustrated in Fig. [5] extracted from [29J. Here the red 
dot-dashed line corresponds to the conventional opacity in the absence of ALPs and the solid 
black line is the average prediction with ALPs. It appears that the average transparency is 
indeed higher than the conventional case at high energy. However, the associated uncertainty 
on the prediction, in other words the variance related to the randomness of the magnetic field is 
such that the envelope includes the conventional case. Because of that fact, if observed without 
ambiguity in the future, such an effect might be seen as an indication for ALP detection but 
could hardly serve as a firm argument for discovery. 
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Figure 5. Transmission of photons with and without ALPs, the ALP case is ploted with 
the envelope corresponding to the variance on the prediction of the transparency effect (figure 
from [29] , with permission of the authors). 



Another limitation comes from the use of a very optimistic value for the extragalactic 
magnetic field. For the ALP effect to significantly affect the opacity, magnetic fields of nG 
strength with Mpc coherence length have to be present in the intergalactic medium. It is actually 
possible to generate such magnetic fields from inflation or QCD phase transition for instance, 
but the required strength is very close to current upper limits. This is illustrated in Fig. [6] 



(from [30]), where the different observational constraints on large scale magnetic field appear 
together with predictions from models (orange thin lines). There the red cross corresponds to 
the typical parameters used in the ALP analyses. It lies in a region that can be seen as fine- 
tuned given the size of the still open parameter space. At the moment it seems invoking such a 
strong magnetic field would be acceptable if the tension in the TeV observations was stronger. 
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Figure 6. Constraints on large 
scale magnetic fields and predic- 
tion of the models. The red 
cross corresponds to the param- 
eters used for the ALP solution 
to the transparency "hint" (figure 
adapted from [30J, with permis- 
sion of the authors). 



A clever way to avoid using intergalactic magnetic fields is to remark that if the source is 
magnetized or embedded in a cluster, then the mixing to ALPs could occur essentially around the 
source. Then, the magnetic field of the Milky- Way can serve as a target magnetic field to convert 
back the ALPs into photons. In that case as well, a strong boost can be expected at high energy, 
as first proposed in [31] and then again in [32j. This effect is used in [33J to estimate possible 
lower limits on the g ia coupling, but again assuming the tension is real between observations 
and models for the transparency. Finally it appears that the transparency observable could be 
used in the future for indication or discovery, if a clear tension was observed. To do so, one 
might wait for the next generation of gamma-ray telescopes such as CTA to have a significantly 
larger sample of sources. 

The problem of having only a few sources can be circumvented by using an energy band for 
which detections are numerous. This has been proposed in [34], where the authors remark that if 
the strong mixing regime is realized, the statistical properties of the observed fluxes from X-ray 
sources could display features distinctive of ALP effects. In that case, sources would be seen 
with fluxes reduced by a factor of 1/3 on average. Of course having no access to the absolute 
intrinsic fluxes, this overall factor is not observable. However because of the random nature of 
the mixing process in astrophysical magnetic fields, the first and second momentum distribution 
should have different shapes compared to the conventional case. In [34] the authors claim the 
observation of anomalous features in the momenta distributions. That result has however shown 
to be questionable in [35] where the effect is claimed to be caused by outliers. Because in that 
case the detection would rely on shapes of distributions, it is difficult to infer a constraint without 
a deeper analysis and this effect is again used to propose a hint. Nevertheless, it illustrates one 
possible use of the stochastic nature of the mixing in astrophysical environments, which is no 
more a limitation but becomes a tool for identifying possible ALP effects. In the following, it 
is shown that a careful study of this randomness can lead to observable effects that are used to 
set constraints on the ALP parameters. 



3. Constraints on ALP parameters from observations of the high-energy sky 

3.1. Effect of the magnetic turbulence 

One peculiar effect of photon/ ALP mixing is the fact that the magnetic field turbulence can 
directly imprint features in the energy spectra from high-energy sources. The exact spectral 
shape one gets at the end is unpredictable, but as shown in [36J the statistical properties of 
the induced irregularities are a prediction of the ALP model. The authors of [37J and [29J 
already noticed that in principle the observed spectra should be very irregular in case of strong 
photon/ ALP mixing, without considering the use of the irregularity as an observable. 

To account quantitatively for the irregularity, the 2 polarizations of the photon must be 
considered, so that the evolution of the system after n domains is given by 

M = II(V exp[-i(E + Mt)s k ) P k ) |^o) , with Ml = P k M k P k 1 (5) 

k 

and 



M k = 
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(6) 



k stands for the fc th domain, of size s k , P k is the rotation matrix between the interaction 
eigenstates and the propagation eigenstates and the matrix M k describes the mixing. The 
indexes k are there to recall that from one magnetic domain to the next, the corresponding 
parameters change due to the different orientations of the magnetic field (B t is the projection of 
the magnetic field on the polarization plane and <j) is the angle that projection makes with one of 
the two photon polarization). m 7 = 4:7van e /m e is the effective mass of the photon propagating 
in a plasma with electron density n e . Examples of spectral oscillation patterns in one domain are 
given in Fig. [7|for different values of S = g ia B t s/2, s being the size of the coherent domain. When 
several domains are considered, the spectrum ends up being very irregular as shown in Fig. [8] in 
the case of an unpolarized beam. For that example, an extragalactic source is considered and 
the magnetic field is typical of that of a galaxy cluster. The top panel of Fig. [8] is the raw signal 
and the bottom panel is the same signal smoothed by the energy resolution of HESS (^15%). 
In that case the critical energy is of order 1 TeV and the effective photon mass is negligible. 

Whereas in the case of the extragalactic magnetic field the naive description of the turbulence 
might be sufficient (essentially because its properties are very poorly known), galaxy cluster 
magnetic fields may deserve a better treatment. The magnetic field in that case is modeled by a 
gaussian field with zero mean and a distribution of modes that is described by a Kolmogorov-like 
spectrum as in Eq. [7| 

k 2 

(SB) 2 oca 2 - -—- . (7) 

The corresponding power spectrum is modeled by a function resembling that of Fig.|9j In galaxy 
clusters, the typical coherence length of the magnetic field is 10 kpc and the strength of the field 
is 1 to 10 fiG. 

3.2. Examples of constraints 

The first example of constraints is from the HESS analysis of PKS 2155-304 [38, 39J, which is 
an AGN located at z — 0.116. For that source, both the extragalactic magnetic field and the 
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Figure 7. Spectral oscillation patterns in 
domains with coherent magnetic field and 
different ALP parameters (figure from [36J). 
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Figure 8. Example of ALP induced 
irregularity in the TeV range (top panel: Raw 
signal, bottom panel: Signal smeared with 
HESS resolution, figure from 
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Figure 9. Typical power spectrum used for 
modeling the magnetic turbulence in galaxy 
clusters. 



cluster magnetic field can be considered. In the first case, as previously, one has to assume 
optimistic values of the magnetic field strength for the irregularity signal to be significant. A 
galaxy cluster is observed around the source, but no magnetic field measurements are available. 
So in the case of the galaxy cluster magnetic field, conservative values for the strength and the 
coherence length are assumed (1 /iG and 10 kpc respectively). As HESS observation ranges from 
hundreds of GeV to a few TeV, from the expression of the critical energy E c it is straightforward 
to see that the typical ALP masses that are probed are of the order of 10 -8 eV. In [38l 139] , 
it is shown that the observed energy spectrum does not exhibit strong irregularities. Then an 
estimator of the irregularity is proposed and numerical simulations are used to exclude sets of 
parameters that lead to significantly too strong irregular behavior. This exclusion has to be 
done on a statistical basis as each realization of the magnetic field turbulence is different. The 



results of the analysis are presented in Fig. [10| The method allows to improve the CAST limits 
in a limited energy range around 20 neV. 

Another possibility is to use a source that lies at the center of a well studied galaxy cluster. 
In that case, the magnetic field properties are derived observationally. This is done by studying 
the Faraday rotation maps of the polarized radio emission from the cluster (see [40J for a review). 
These studies allow in principle a determination of the full turbulence power spectrum, yielding 
the intensity of the magnetic field, its coherence scale and the slope of the turbulence spectrum. 
A very well studied cluster is Hydra, for which a strong X-ray source is present at the center 



(Hydra A) [H] . In [12], X-ray data from the Chandra satellite are analyzed in order to derive 
constraints on ALP parameters. In the case of X-rays from Hydra, the diagonal terms in the 
matrix of Eq. [6] can be simplified. Indeed the pair production related opacity is irrelevant in the 
case of X-rays (so r = 0), and the trace of the matrix is dominated by the effective photon mass 
for m a < 10 -11 eV. So the constraints are expected to extend to arbitrarily low ALP masses 
below that value. In [42J, the irregularity is estimated by performing x 2 tests when deriving the 
energy spectrum with a forward folding method. ALP parameters yielding a too high level of 
irregularity compared to the data are excluded. The corresponding exclusion curve is displayed 
in Fig. [TT] It turns out this analysis improves the previous constraints in that mass range from 
the non-observation of gamma-rays associated with SN 1987 A [^3] . 
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Figure 10. HESS exclusion contours from Figure 11. Exclusion contours from the 
the analysis of PKS 2155-304 energy spectrum analysis of X-ray data from Hydra A (Figure 
(Figure from [38l|39j). from 02]). 



4. Outlook 



The study of the high-energy Universe is potentially a nice way to search for axion-like particles. 
The problem of the transparency of the Universe to gamma-rays can provide a interesting 
observable, that requires nevertheless the observation of a large number of TeV sources to be 
robust. This can be achieved with the next generation of Cherenkov telescopes such as CTA. 
The photon/ ALP mixing in astrophysical sources is intrinsically a stochastic process. That 
fact makes difficult the use of the transparency observations to derive constraints on the ALP 
parameters. It is noted however that the turbulence of the astrophysical magnetic fields has the 
effect of inducing irregularities in the energy spectra of sources. The statistical properties of the 
induced irregularity can be predicted and are used to set limits on the ALP coupling to photons. 
Because the method is insensitive to the polarization, these constraints go beyond classic ALPs 
and apply to both FF and F 2 types of couplings. Two examples of limits are given in the case 
of TeV and X-ray observations of high-energy emitting sources inside clusters of galaxies. 
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